Decays of spacelike neutrinos* 
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In this paper we consider the hypothesis that neutrinos are fermionic tachyons with helicity 1/2. We 
propose an effective model of interactions of these neutrinos and analyze dominant effects under the 
above hypothesis: the three-body neutrino decay, the radiative decay and briefly discuss /3-decay 
with tachyonic antineutrino. We calculate the corresponding amplitudes and the mean life time 
for decays, as well as the differential decay width. In addition we apply the neutrino three-body 
decay to the solar neutrino problem and we obtain a remarkable change of energy spectrum of the 
neutrino flux at the Earth ground (assuming no neutrino oscillation). We compare the results with 
the Standard Solar Model predictions and solar neutrino experiments to estimate an upper bound 
for the taonic neutrino. 



I. INTRODUCTION 

The indication for negative square of electron and 
muon neutrinos four- moment a in all the recent experi- 
ments ji]]^] requires a theoretical or methodological ex- 
planation. However, up to now, there is no satisfactory 
way to understand this phenomenon. The most pop- 
ular solution of this problem is a possible presence of 
molecular effects in the final state; however they are not 
understood well j^]. Another proposed possibility is an 
existence of a hidden anomalous long range interaction of 
neutrinos [0 or a presence of a repulsive potential origi- 
nating in coherent interaction of neutrino with the daugh- 
ter atom ||. Also some methodological misinterpreta- 
tions can cause the problem, as it was suggested by Roos 
and Khalfin || ; see also (7| for a discussion of the negative 
mf, puzzle. However, if the negative squared antineutrino 
mass in tritium beta decay is genuine, then it is hard to 
understand this phenomenon using conventional particle 
physics ideas. Moreover an anomalous electron/muon 
composition of flavor of atmospheric neutrino flux |}]|^] 
as well as the solar neutrino problem is a puzzle, which 
is not solved to the end, even if we take into account 
possibility of flavor oscillations and Mikheyev-Smirnov- 
Wolfenstein (MSW) mechanism §-|l|. 

In this paper we continue investigation of an unconven- 
tional possibility that neutrinos are fermionic tachyons. 
This hypothesis was proposed firstly by Chodos et al. 

In the papers by Rembielinski |HJ a causal clas- 
sical and quantum theory of tachyons was elaborated 
and next applied by Ciborowski and Rembielinski to ex- 
planation of the results of tritium /3-decay experiments 
p5|^6| . In the paper [|l5| physical consequences of the 
three body decay vg — > vivyvgi was also discussed. This 
paper is devoted to systematical calculations of decay 
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rates for the two processes with participation of neutri- 
nos: the radiative decay — > v^y and the three body 
decay vg — > vgDgiVgj, both conserving the lepton flavor; 
we briefly discuss also the /?-decay. Notice, that the first 
two reactions are kinematically admissible for tachyonic 
neutrinos only. It is remarkable that, as was stressed in 
|]l5| , the emission by neutrino vg a vt'Vt' pair can be an 
additional, qualitatively similar pro cess to flavor oscilla- 
tions. As we will see (Sec. |v| and VH), the calculated 
mean times of life for the aforementioned decays take 
reasonable values in the regions of experimentally avail- 
able energies. Notice however that applicability of these 
results to the description of the time evolution of a neu- 
trino flux is restricted at this level because both of the 
decays form cascades. Therefore they demand further in- 
vestigation as nontrivial stochastic processes leading to 
different than Geiger-Nutall decay law. Moreover, the 
mostly hopefull program is to involve both oscillations 
and three body decay. The corresponding calculations 
for solar and atmospheric neutrinos as well as for the 
supernova SN1987A neutrino flux are now in progress. 
However we give here an estimation for the upper bound 
of the tau neutrino mass (if it has tachyonic character) 
by means of the solar neutrino data. 



II. PRELIMINARIES 

In this chapter we give a brief presentation of the causal 
theory of tachyons proposed in the paper [S3 . 

It is a common opinion that existence of space-like par- 
ticles is in contradiction with special relativity. This is 
because of the Einstein-Poincare (EP) causality viola- 
tion and consequently because of a number of very se- 
rious difficulties like causal paradoxes, inconsistent kine- 
matics with impossibility of a covariant formulation of 
the Cauchy problem, unbounded from below energy spec- 
trum and many others. Even more problems arise on the 
quantum level, like the problems with construction of co- 
variant asymptotic spaces of states, vacuum instability, 
etc. 

However, as was shown in the paper E3| it is possi- 
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ble to agree special relativity with the tachyon concept. 
The main idea is based on the well known fact that the 
definition of the coordinate time depends on the synchro- 
nization scheme [|T^[l9| which, in turn, is a convention 
related to the assumed one-way light velocity (only the 
average of the light velocity over closed paths has an 
operational meaning; in special relativity this average is 
frame independent). Taking into account this freedom, 
it is possible to realize Poincare group transformations 
in such a way that the constant-time hyperplane is a 
covariant notion In that synchronization (named 

Chang Tanghcrlini synchronization scheme and denoted 
as CT synchronization — see |l4j and references therein) 
it is possible to overcome all the difficulties which appear 
in the standard approach. 

Let us stress main features of this nonstandard scheme. 
Firstly, it is fully equivalent to the standard formula- 
tion of the special relativity if we restrict ourselves to 
the timelike and lightlike trajectories, i.e. if we exclude 
tachyons. If we admit space-like trajectories (tachyons 
are included) then one inertial frame is distinguished as 
a preferred frame^J. Thus the relativity principle is bro- 
ken, however the Poincare symmetry is still preserved in 
that case. The proper framework to this construction is 
the bundle of Lorentzian frames; the base space is sim- 
ply the space of velocities of these frames with respect 
to the preferred frame. For this reason the transforma- 
tion law for coordinates involves velocity of distinguished 
frame. The preferred frame can be locally identified with 
the comoving frame in the expanding universe (cosmic 
background radiation frame), i.e. the reference frame of 
the privileged observers to whom the universe appears 
isotropic p0| . 

Now, the velocity of the Solar System deduced from the 
dipole anisotropy of the background radiation is about 
350 km/sec so it is almost at rest relatively to the 
preferred frame. Therefore, with a good approximation, 
we can perform calculations in the preferred frame. 

To be concrete, the Lorentz group transformations in 
the mentioned bundle of frames have the following form 



x' = D(A,u)x , 
v! = D(A,u)u, 



(1) 



where for rotations D(R,u) has the standard form while 
for boosts it reads 



D(W,u) 



-w 



W&W 1 



W <g> U T U° 



(2) 



Here W 11 is the four- velocity of (x 1 ) frame seen by an 
observer in the frame (x) while u M is the four-velocity of 
the privileged frame as seen from the frame (x). Notice 
that the time coordinate is rescalled by a positive factor 
only. The transformations (^) leave invariant the metric 
form 



ds 2 = g^ v (u)dx^ 'dx v 



with 



9{u) 



-fT 

u u 



-I + u<S>u T (u°) 2 



(3) 



(4) 



Interrelation with coordinates in the Einstein-Poincare 
synchronization (x£) is given by 



x° ± u°u ■ x . 



x E 



(5) 



However the corresponding interrelations between veloc- 
ities and v obtained from (|^) are singular for superlu- 
minal velocities. 

Now, on the quantum level, the corresponding quan- 
tum mechanics and quantum field theory can be formu- 
lated with help of the bundle of Hilbert spaces associ- 
ated with the aforementioned bundle of frames. The 
Fock construction can be done in this framework and 
unitary orbits can be classified For usual particles 
they coincide with the standard unitary representations 
of the Poincare group. For tachyons they are induced 
from SO{2) group (instead of 5*0(2, 1)) and labelled by 
the helicity. The following two facts, true only in CT syn- 
chronization, are extremely important for quantization of 
tachyons: 

• Invariance of the sign of the time component of the 
space-like four-momentum, i.e. e(k°) — inv, which 
follows from the transformation law (|l|)-(||). 

• Existence of a covariant lower energy bound; in 
terms of the contravariant space-like four-momen- 
tum k* 1 (k 2 < 0) this lower bound is exactly zero, 
i.e. k° > as in the timelike and lightlike case. 

This is the reason why an invariant Fock construction 
with a stable vacuum can be done in this case on the 



quantum level |14|. Notice that the measure 



dn(k, k) = d 4 k 0(k°) 5{k 2 + k 2 ) 



(6) 



is Poincare covariant in the CT synchronization scheme. 



III. FERMIONIC TACHYONS WITH HELICITY 

A = ±1/2 



*We would stress that the notion of preferred frame used 
through this work is not the one used in the search of breaking 
the Lorentz symmetry. 



To construct tachyonic field theory describing field ex- 
citations with the helicity ±1/2, we assume that our field 
transforms under Poincare group like bispinor (for dis- 
cussion of transformation rules for local fields in the CT 
synchronization see |ljj2^]); namely 
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U(A)iP(x, u)U(A~ 1 ) = SiA- 1 )^',^) , 



(7) 



where S(A) belongs to the representation D^ D°5 of 
the Lorentz group. Because we are working in the CT 
synchronization, it is convenient to introduce an appro- 
priate (CT-covariant) base in the algebra of Dirac matri- 
ces as 



(8) 



where 7^ are standard 7-matrices, while T(u) is deter- 
mined by the Eq. (||). Therefore 



{r,l u } = 2g^{u)I. 



(9) 



However, notice that the Dirac conjugate bispinor ip = 
^7g- Furthermore 7 s = — i£my CT A7^7 t/ 7 cr 7 A /4! = 7|. 
Let us consider the Dirac-like equation)] [ fL4[ 

(7 5 (i 7 <9)-K)7/; = 0. (10) 

which can be derived from the Lagrangian density 

£ = 1 /i( 7 5 (i 7 a)- K )^. (11) 

Dirac equation (|l0|) imply the Klein-Gordon equation 

(g^(u)d^- K 2 )ij = 0, (12) 

related to the space-like dispersion relation k 2 = —n 2 . 
The equation (JTTJ) is analogous to the one by Chodos et 
al. |p^| Dirac-like equation for tachyonic fermion. How- 
ever, contrary to the standard EP approach, it can be 
consistently quantized in the CT scheme if it is supple- 
mented by the covariant helicity condition |14| 



X(u)tp(u, k) — \ip(u, k) 



with A given by 



X(u) 



y/(Pu) 2 ~ P 2 1 



(13) 



(14) 



where 



W» = 1/2 e^ AT J aX Pr 

is the Pauli-Lubanski four-vector. This condition is quite 
analogous to the condition for the left (right) bispinor in 
the Weyl's theory of the massless field. It implies that 
particles described by ^ have helicity —A, while antipar- 
ticles have helicity A. For the obvious reason in the fol- 
lowing we will concentrate on the case A = 1/2. 

Notice that the pair of equations ([iTf), ( fl3| ) is not in- 
variant under the P or C inversions separately. 



Now, in the bispinor realization the helicity operator 
A has the following explicit form |i~4j 



A(«) = 



7 5 [i7<9, zt 7 ] 



(15) 



— 1/2 

where the integral operator ((iud) 2 + □) in the co- 
ordinate representation is given by the well behaving dis- 
tribution 



= (2^- 4 / d%e(u V yv x ((up) 2 - P 2 y 1/2 



(16) 



The free field ip(x,u) has the following Fourier decom- 
position 

ip(x, u) = (2tt)- 3 / 2 J d 4 k 6(k 2 + K 2 )9(k°) 

x [w(k,u)e ikx tf(k) +v(k,u)e- lkx a(k)] . (17) 

The creation and annihilation operators of a tachyonic 
fermion (a) with helicity —1/2 and an antifermion (b) 
with helicity 1/2 satisfy almost standard canonical anti- 
commutation relations; the nonzero ones are 



[o(fc),at(p)] + =2^S(k-0, 
[b(k),tf(p)] + =2u % 5(%-2). 



(18) 
(19) 



Here ui^ — k° > is the positive root of the dispersion 
relation k 2 — —k 2 and k denotes covariant components 
fej, i = 1, 2, 3 of k. The amplitudes w and v satisfy 

(K + 7 5 fc 7 ) w(fc,u) = 0, (20) 

(l+ 75 j fc7 '" 7] )w(k,u)=0, (21) 

(K~-/ 5 kj)v(k,u) = 0, (22) 

fl + ^4) «(*,«) =0. (23) 

Furthermore, the projectors ww and vv read 



^Hereafter wy = u M 7 M , ud = u^d^, 7 9 = j^d^ 



W (k,u M k,u) = ^l[l-^^), (24) 

v(k,u)- V (k,u)=- K -±pl(l-^^Y (25) 
2 y 2yJq 2 + K 2 ) 

The above amplitudes fulfill the covariant normalization 
conditions 

w(k, u)j 5 u^fw(k, u) = v(k, u)^ 5 wyv(k, u) = 2uk , (26) 
u)(fc ,r ,M)7 5 u7w(fc,u) = 0, (27) 

where the superscript tt denotes the space inversion of k. 
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It is easy to see that in the massless limit k 
Eqs. (|20l)-(f23|) give the Weyl equations 



the 



kjw — kjv = . 



7 5 w 



-w . 



Now, the normalization conditions (p6|) -(p7 [) together 
with the canonical commutation relations (|lq), (|l9|) guar- 
antee the proper work of the canonical formalism. In 
particular, starting from the Lagrangian density ([ll]) we 
can derive the translation generators 



P„ 



d 3 k{2uj j :y 1 k tl (a t (fc)a(fc) + {k)b{k)) . (28) 



Thus we have defined a consistent Poincare covariant free 
field theory for a fermionic tachyon with hclicity —1/2; 
the Weyl's theory for a left spinor is obtained as the k — > 
limit. For more details of this construction see fl4|| . 



IV. DYNAMICS 

To formulate a dynamical model with tachyonic neu- 
trino we can use, on the tree level, the phenomenologi- 
cal form of the weak interaction Lagrangian, namely the 
current-current Fermi Lagrangian. The only condition 
we take into account is that the massless limit n — ► of 
our model should coincide with the standard four-fermion 
interaction. This leads to three natural possibilities for 
the part of the lepton charged current containing tachy- 
onic neutrino: 



I(x) 7 »±(l--y 5 )v(x) 
i_{x)^v{x) + h.c. 
£(x)j 5 ^v(x) + h.c. 



h.c. — chiral coupling, 
— helicity coupling, 
— 7 5 coupling, 



and similarly for the neutral currents. 

Recall that v(x) satisfy the helicity condition (Jl^) with 
A = 1/2 so in the K — > limit v — > v^. 

However, if we consider higher order processes (like 
v — > wy) it is necessary to have a more sophisticated 
model like the Salam- Weinberg one. Let us notice firstly, 
that the Lagrangian for the massless Dirac field, under 
assumption of CP conservation, can be written in two 
different forms ipijdip or tp^i'jdip. Both forms lead to 
the same massless Dirac equation ijdip = because 
of invertibility of 7 s matrix. Therefore generation of 
a mass via Yukawa coupling with Goldstone fields can 
lead to massive fermions or fermionic tachyons, because 
of creation in the Lagrangian of Dirac -0 (i-fd — m) i/j or 
■0 (j 5 ijd — to) V> term respectively. Thus, because we 
have decided to treat neutrinos as tachyons, we start with 
the following free massless Lagrangian for a fixed lepton 
generation (y, £) 



(29) 



After splitting v and I into left-handed and right- 
handed pairs vl,r, £l,r, it is easy to see that the (com- 
pact) invariance group of the Lagrangian Cq is exactly the 



(SU(2) x U(1)) L x U(1) R group! Therefore, taking into 
account the electric charges of v and £, we conclude that 
the symmetry group of £ is lepton number X (U(l)y x 
SU(2)ij and the left-handed fields Vl and II must form 
the weak dublet while the right-handed vr and £r are sin- 
glets of the weak group. Therefore Cq leads uniquely to 
the standard Salam- Weinberg choice of the weak symme- 
try group and their realization on the lepton fields (recall 
also that in the k — > limit v(n) — * u(0) — vl(0), i.e. 
vr(k) — > ^r(O) = 0). Now, under the standard choice 
of the Goldstone fields (as the U(1) Y x 5C/(2)/ dublet) 
as well as the Higgs potential and Yukawa coupling (for 
neutrinos like for neutral quarks) and after gauging the 
weak group we obtain the final Lagrangian 



C ^fermions + int. 

~t~ (-^gaugc bosons ~t~ ^Higgs + int. 



(30) 



The part of C in the parentheses has the standard form 
(we choose convention as in [p3| ) while the £f C rmions + int. 
reads 

Acrmions + int. = V^i^dv - KVV 

+iijd£ - mil + eA^H 



f) 



2 cos 9w 



Zyfrf (9V - 9A1 5 ) £ 



\pl 2 \/2 2 

? J v ((to - k)I + (to + k) 7 5 ) IG+ 

2\l2mw 

£([m - k)I - (to + k)7 5 ) vG~ 



2V2: 

ig_ 

9 

2mw 



raw 



{kvv + mli) H . 



(31) 



Here gy = 1/2 — 2 sin 2 6w, 9 'A = 1/2, g — weak coupling 
constant, 9w — Weinberg angle, G^ 1 and G° — Goldstone 
bosons and H — Higgs field. The corresponding Feynman 
rules are listed in the Appendix 

The following remarks are in order. First of all, the 
derivation of the Lagrangian ([5(j|) with £f er mions + int. 
given by ([H]) is rather heuristic one; the point is that the 
right component of v, ur goes to zero with k — > 0. This 
means that the mass generation mechanism for tachyons 
is not well defined by Yukawa coupling if we start with 
the massless Lagrangian. Intuitively, it seems that the 
mechanism of mass generation for the tachyonic neutri- 
nos should be caused rather by the interaction with grav- 
itational field. This is related to our natural assumption 
about coincidence of the tachyon preferred frame with the 
comoving frame. The next problem is the hclicity con- 
dition (O) which is necessary for consistency of the free 
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FIG. 2. Effective four-fermion vertex in neutrino three- 
-body decay vt — > vgv^vgi 
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FIG. 3. The mean time life r«< in the decay vg — > vgvgivgi 
as the function of the pair neutrino (i^/) mass for three 
different values of the mass k of the decaying neutrino vg . The 
plot is given for the three values of the total energy E — k° 
in the preferred frame 



tachyon theory. Of course it can be gauged by introduc- 
tion of the covariant derivatives and causes additional re- 
lations between multipoint Green functions of the theory 
In such a case perturbation consistency and renormaliz- 
ability of this model should be proved. We shift these 
unanswered questions to further investigations and here 
we will treat the Lagrangian (^0|) as an effective reason- 
able approximation of a more realistic theory. Because in 
the following we do not analyze processes with intrinsic 
neutrino lines, this approximation seems to be correct. 

Notice that the Feynman rules derived from (|30|)-(|3l|) 
(Appendix [X]) , lead to the effective chiral coupling in the 
Fermi Lagrangian. Therefore, in the following all the cal- 
culations on the tree level will be done for this coupling. 
However, we will comment results obtained with help of 
the helicity and 7 5 couplings too. 

Tachyonic neutrino is in general unstable even if the 
leptonic flavors are conserved; namely decays of the form 
tachyon — > tachyon + sth. are kinematically admissible. 
If we take into account the experimental upper bounds 
on the neutrinos masses jj| we can select two dominant 
processes of this type, namely: three body decay vg — > 
Vn + vgi + &£/ and radiative decay vg — ► v g + 7. 

V. NEUTRINO DECAY vg -> v,vg,vg, 

This process can be analyzed on the tree level. The 
corresponding amplitude is the sum of the graphs given 
in Fig. ^ and leads to the effective four-fermion vertex in 
Fig. g. Consequently 

M = V2G F v e (p)r^v e (k)v e ,(l)T^w e ,{r) , (32) 

where the momenta k, p, I and r correspond to the incom- 
ing neutrino vg and to the outgoing neutrinos vi, vy and 
vgi respectively. The masses of v and v# are denoted by 



k and [i respectively. In our case the Feynman rules (Ap- 
pendix^) lead to the chiral coupling T M = (1 + 7^)7^/2. 
Thus, by means of the polarization relations (E4j)-(Eq), 
we obtain in the preferred frame 

or 2 

\M\ 2 = -^-(\k\ + k )(\p\ + p°)(\l\ + l ) 
\k\\p\\l\\f\ 

x(\r\+r°) (\p\\l\-P-l) (\kM-k-f) (33) 
where from dispersion relations \k\ — k 2 + (fc ) 2 , |p| = 

V« 2 + Cp ) 2 , \i\ = vV + (i ) 2 , M = vV + (r ) 2 . 

The general forms of the square of the amplitude for all 
couplings are given in the Appendix |b|. The qualitative 
behavior of the decay rate 

r«, = ((2^) 5 2fc°) _i y d A P d A id A r \m\ 2 e(p )e(i )e(r ) 

x S(p 2 + k 2 )5{1 2 + n 2 )S(r 2 + n 2 )5 4 (k —p — l — r) (34) 

is rather similar for all three couplings. Here we present 
in Fig. [| the numerical calculations for the mean life time 
in the case of chirality coupling (according to our gauge 
model) done in the preferred frame [u = (1,0)). For 
the helicity coupling the corresponding results were given 
in p5| . The mean life time of the tachyonic neutrino 
in this process is determined by its energy E — k , as 
measured in the preferred frame, and masses of all three 
neutrino species. As we see from the Fig. || the life time 
tu' {£ 7^ £') corresponding to the partial width Ta', is 
a decreasing function of the initial energy and mass of 
the neutrino vy. For illustration we present in Fig. [| the 
energy spectrum of the final neutrino vg (after decay) 
assuming its mass as k = 10 eV and for the vy , taking 
(i = 0.15 MeV and fi = 1.125 MeV. The case 1 = 1' 
differs by the corresponding combinatorial factor. It is 
important to stress that: 
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FIG. 4. Differential decay rate dF/d£ for the three-body decay u e — > v e v^ii^ for v e mass re = 10 eV as the function of the 
outgoing neutrino v e energy £ = p° for £' = fi (on the left), r (on the right), with /i M = 0.15 MeV, /i T = 1.125 MeV respectively; 
the incoming neutrino energy k° = 10 MeV 



• The decay vi — > v^vgiugi has character of an emis- 
sion by a neutrino-antineutrino pair P^/ ; as 
a consequence the resulting neutrino flux contains 
neutrinos and antineutrinos of all possible flavors. 
Therefore this process can simulate the neutrino 
flavor oscillations. Notice that for massive and 
massless neutrinos such a decay is kinematically 
forbidden. 

• The above process is repeated and forms a cas- 
cade. Therefore, as was mentioned in the Intro- 
duction, it demands more sophisticated treatment 
as a nontrivial stochastic process. However after 
the neutrino energy degradation it is slowing down 
with subsequent decays; notice that, for low ener- 
gies in the preferred frame, neutrinos are almost 
stable (Fig. H). 



VI. THREE BODY DECAY AND THE SOLAR 
NEUTRINO PROBLEM 

As is well known, predictions of the Standard Solar 
Model (SSM) are in a strong disa gree ment with all oper- 
ating solar neutrino experiments p4]-p6[ . 

Firstly, the chlorine radiochemical experiment (thresh- 
old ~ 0.81 MeV) gives measured event rate which is a fac- 
tor 3.6 less than the theoretically predicted by using SSM 
with helium diffusion^]. The neutrino flux detected in the 
chlorine experiment is from the high energy 8 B neutrinos 
with a significant contribution of 7 Be neutrinos; contri- 
butions of pep and CNO neutrinos are not so significant, 



*SSM with helium diffusion are mostly favored by the helio- 
seismological experiments p7| . 



p8j . Secondly, the water experiment (Kamiokande; ap- 
paratus threshold ~ 7 MeV) detects 8 B neutrinos and 
gives about 0.4 part of the predicted 8 B neutrino flux 
|||. Finally, the gallium experiments (GALLEX, SAGE) 
detects mostly basic pp neutrino flux as well as pep neu- 
trinos, but also 8 B, 7 Be and CNO fluxes contribute to 
the measured event rate p8| . Observed event rate in this 
experiment is about 0.52 of the SSM prediction pfjj . 

Let us check applicability of the tachyonic neutrino 
three-body decay to the solar neutrino problem. Al- 
though vacuum oscillation or MSW mechanism are possi- 
ble in this case too, we assume here that neutrino deficit 
is caused by the neutrino decay only. As we will see, 
this assumption allows us to estimate the taonic neu- 
trino mass (if it has tachyonic character) as less than 
1.125 MeV. In the following we use the neutrino fluxes 
as predicted in papers by Bahcall and Pinsonneault |28j . 

As it is evident from the Fig. in the electron neu- 
trino (k c ~ 10 eV) decay, for energies 1-10 MeV and 
time Sun-Earth distance/c ~ 500 sec, the mass of the 
resulting vy should be of the order 1 MeV. This excludes 



decays v e 



and v e 



as possible sources 



of the investigated effects (recall the Jackelmann A and 
B solutions |2^] for the muon neutrino mass ^0.37 MeV 
and ~ 0.15 MeV respectively). We show the dependence 
of the mean life time of v e on energy for these cases on 
Fig. |^. Therefore in the following we will consider the 
contribution of the v e — > v e v T v T decay only. Our aim is 
to analyze the cascade process 



V e — > v e v T v T 
i 

V e V T V T 

I 

Ve V T V T 

I 

taking into account that probabilities of the possible sub- 
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FIG. 5. The mean time life r e i of electron neutrino in 
three-body decay into electron neutrino and neutrino-anti- 
neutrino pair of type I for t = e (mass k = fi e = 10 eV), i — ft 
(A) and (B) (masses ^ M = 0.37 MeV and ^ M = 0.15 MeV cor- 
responding to the Jackelmann A and B solutions respectively 
§9|), and for i = r (mass fi T = 1.125 MeV) 

sequent decays v T — ► v T VeV e and v T — > v T VeP e are not 
significant in the investigated energy regime. 

This greatly simplifies treatment of this stochastic pro- 
cess leading to the following evolution equation for the 
electron neutrino flux: 

n ^dP~ - ~™ E) + £ dE ' d -^r ME ' ] (35) 

Hereafter $t(E) denotes the electron neutrino flux of en- 
ergy E after time t, T E is the total decay rate for the 
initial neutrino of energy E and dT e> / dE is the differen- 
tial decay rate of the decaying neutrino v e with energy 
E' into outgoing electron neutrino with energy E. The 
equation ( p5[ ) has the solution as a formal power series 

$ t = e tF / R <fo (36) 

with F(E, E') = -T E 5(E - E') + 9(E' - E) dT E ,/dE. 

The numerical estimation of the spectrum of decay- 
ing neutrinos with respect to the heaviest neutrino mass 
(y T ) gives the best fimfor fi T = 1.125 MeV. In this case 
the neutrino spectra p^p0| -p2[ lead to the results showed 
in the Table ffl. It is evident that, under assumption of 
the correctness of the SSM predictions, the three-body 
neutrino decay does not solve the problem to the end. 
However we can see that the problems mentioned above 
can be reduced to the overflow of 7 Be neutrinos only. In- 
deed, if we set 7 Be flux zero, we obtain approximately 
$ pp ~ $pp M and $8 B — 0.4$f| M , consistently with ob- 
served values. We would point out, that in (331 it was 




FIG. 6. Normalized (to unity) energy spectrum of S B neu- 
trino flux registered on the Earth from SSM predictions and 
our model; the masses are assumed to be fi e = 10 eV, 
fir = 1.125 MeV 
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FIG. 7. Energy spectrum of 7 Be neutrino flux from SSM 
predictions and our model; the masses are assumed to be 
fi e = 10 eV, fir = 1.125 MeV 

argued that a slow mixing of 3 He down into the solar core 
can cause significant reduction of the 7 Be flux. On the 
Fig. | and Fig. | we compare calculated by ourselves and 
predicted by SSM 8 B and 7 Be neutrino fluxes. In princi- 
ple, the change of the shape of 8 B neutrino flux (if it actu- 
ally changes) , can be measured in the SuperKamiokande 
experiment. 

Finally, let us stress that the above results was ob- 
tained under the assumption of the absence of possible 
neutrino oscillations But they can be used to an estima- 
tion of the upper bound for the taonic neutrino mass. 
We intend to return to the more general situation (with 
oscillations) in the forthcoming paper. 



s For the helicity coupling we obtain for fi T the value 
0.95 MeV. 
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TABLE I. Solar neutrino fluxes — comparison of our model and SSM 





2 H 




37 C1 






71 Ga 






our model/SSM 


our model 


SSM 


our model/SSM 


our model 


SSM 


our model/SSM 




% 


(SNU) 


(SNU) 


% 


(SNU) 


(SNU) 


% 


pp 










69.30 


69.69 


99 


pep 




0.19 


0.21 


93 


2.71 


2.86 


95 


Be 




1.38 


1.41 


98 


37.37 


37.82 


99 


8 B 


41 


2.83 


7.01 


40 


7.35 


16.16 


45 


13 N 




0.10 


0.11 


94 


3.77 


3.88 


97 


15q 




0.34 


0.36 


95 


6.13 


6.33 


97 


17 F 




0.0041 


0.0043 


95 


0.073 


0.076 


97 


total 


41 


4.84 


9.10 


53 


126.70 


136.82 


93 


our model — 7 Be 


41 


3.46 


9.10 


38 


89.33 


136.82 


65 



VII. RADIATIVE DECAY OF THE TACHYONIC 
NEUTRINO v t -> vi + 7 

This process can play an important role in cosmology; 
in particular, as was noticed in |m] [36|, it can affect the 
cosmic background radiation in the submilimeter region, 
leading to its distortion. 

Let us denote the four-momenta of the initial and fi- 
nal neutrinos by k and p respectively, their masses by 
k, while photon has the four-momentum q = k — p. 
The kinematics of this process is very simple; namely 
k 2 = —K 2 ,p 2 = — k, 2 and for the photon on the mass 
shell q 2 = (k — p) 2 = so kp = —k 2 , qk = qp = 0. Let us 
consider this kinematics in the preferred frame u = (1,0). 

It is immediate to see that cosO = k ■ q/(\k\\q\) = 
r /2 

k° ((k ) 2 + k 2 ) 1 . Therefore the out going photons are 
emitted on the cone determined by the angle 28. Notice 
that for fc°>K, cos 9 w 1 so the process is approximately 
collinear in that case. 

The amplitude has the form 



(37) 



where is the polarization vector for the photon with 
helicity A. The fi M (p, k) should satisfy transversality and 
neutrality conditions. The first one follows from the elec- 
tromagnetic current conservation and has the form 



q^{p, k) = 



(38) 



while neutrality condition means that the neutrino charge 
is zero and consequently the corresponding matrix ele- 
ment vanishes: 



fF(fc,/s) = 0. 



(39) 



By means of the equations of motion for the spinors v 
and v (p2[), the most general form of the f2^(p, k) reads 



W(p,k)=v(p){(F v + F A7 5 ) 1 >* 

+ (G lv + iG 1Al 5 )k» 

+ {G 2V +iG 2Al b )p»}v{k), 



(40) 



where the form factors F and G are in general functions 
of q 2 . Notice that the form factors G{ A control CP non- 
invariant terms in the Eq. (^0|). Now, taking into account 
the transversality ([38|) and neutrality (|39j) we obtain the 
following constraints 



F v = 0, 

F A = k(Giv + G 2 y) . 
(G2A — G\ A )q 2 = , 
(G 2V - G lv )q 2 = . 



(41) 



Therefore, because of the analyticity of Giv and Gt A in 
the point q 2 = 0, we obtain 

n"(p, k) = -hu{p) [Y\ 7 "] q v (G v + iG Al 5 )v(k) , (42) 

where we denote GV = G\y = G 2 v, G A = G\ A = G 2A 
and the equations of motion ( p2] ) was used again. There- 
fore the square of the amplitude Af A , after summation 
over final polarizations of the photon (Appendix ^), 
takes the form: 



|M| 2 -4 K 4 (|G y | 2 + |G A |' 
{uqf 



y/(u(k - q)) 2 + K 2 ^J(uk) 2 + K 2 



(43) 



where we have used the form of the polarization operator 

In the framework of our SU(2) x U(l) gauge model, 
the diagrams contributing to the the process Vt — ► vg + 7 
to one-loop order are shown on Fig. [§. 

Now, the transversality and neutrality conditions (|3S 
( p9| ) can be treated as renormalization conditions. By 
means of the Feynman rules listed in the Appendix A we 
can easily calculate the form factors G A and Gy: 



G A = 0, 



(44) 



as we expected because of the CP invariance of this 
model, while 



Gv - 64^ [Ia 



le + If] 



(45) 
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FIG. 8. Contributions to the process vt — *• vt + 7 to one-loop order 



im 2 w ' 



where (see Appendix |^) , under the experimental condi- 
tions for lepton masses m <C mw and for corresponding 
neutrino masses k <C mw, 



h 
lb 

h 
h 



6m w ' 

2rn 2 (mw 
— t- In 



5m 2 + k 2 
6m^ 



2m 2 w 



Thus in this approximation I c + 1^ can be omitted and 
the final form of Gy reads 



G v = 



eg 2 K 
96ir 2 mw 



kGf I ol 

6tt V 27T 



Consequently 



and 



|M| : 



aK 6 G 2 



18tt 3 ^( M fc) 2 + K?y/(u(k - q)) 2 + k 2 



(46) 



(47) 



(48) 



Notice, that (|47j) implies that the magnetic moment \x v 
of the neutrino v is given by 



kGf I a Gpm e K 



3tt V 2vr 3V2tt 2 



(49) 



where /is is the Bohr magneton and m e the electron 
mass. (N.B. the value of \x v is of the same form as for 
the massive neutrino; the difference lies in the numerical 
factor only [p7|). To calculate the decay rate we integrate 
\M\ over the phase space: 

r=(8f¥)- 1 / d\d A P \M\ 2 e{ q a )6{p Q ) 

xS{q 2 )5{p 2 + K 2 )S(k-p-q). (50) 

We shall do it in the preferred frame, that is we put 
u = (1,0). Inserting \M\ given by ( fig ) we obtain the 
energy spectrum of emitted photons as: 



dT 



aG 2 F K 6 (q Q ) 2 e{k° - q°)0(q°) 



dq° 2887T 4 fc° ((k ) 2 + k 2 ) v /(fc°-9°) 2 + K 2 
Finally, the total decay rate reads: 



(51) 



an 6 G 2 F 



5767T 4 fc° {{k ) 2 + K 2 



4nk° -3fcV(fc ) 2 + ^ 2 



+ (2(fc ) 2 - K 2 ) In V(fc°) 2 + « 2 + fc" 



(52) 
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TABLE II. Upper bounds on the neutrino magnetic mo- 
ment 



from ( 


49 


) from |Q 


Me 

Mm 

i'r 


< 1.4 x lCT 18 ^ 

< 2.4 x lCT 14 /i S 

< 3.4 x 10- 12 /is 


< 1.8 x Kr 1( Vs 

< 7.4 x 10~ 1 Vs 

< 5.4 x 10"Vs 




A 

FIG. 9. Mean life time r in the radiative decay v — > ^7 as 
the function of the neutrino mass ft and the initial energy k°. 
Here £ « 2.6899 



Note that the experimental upper bounds on the neu- 
trino masses applied to the magnetic moment formula 
( [49| ) give results which are in agreement with the exper- 
imental data (cf. Table 0). 

Now, for illustration we present the mean life time r = 
Ti/T for this process as a function of the energy k in 
Fig. |^. From (|5^) it is easy to see that for fc° = with 
£ ~ 2.6899 the mean time of life has a minimum. The 
minimal value r m j n as the function of the mass k reads 

902tt 4 ^ _ 5 

and is presented in Fig. |l^. Finally, in Fig. |ll| differential 
rate dT/dq as the function of the photon energy q° and 
the initial energy A; is given. 



10 40 




10 1V 6 10- 5 10- 4 0.001 0.01 0.1 1 10 100 

k (MeV) 

FIG. 10. Minimal mean time of life T m i n for the radiative 
decay v — + 1/7 as the function of the neutrino mass k 



k" = 0.5k -- 




FIG. 11. Differential rate dT/dq° for the radiative decay 
v — > vy as the function of the photon energy q° close to the 
minimum of the mean time of life 



VIII. /3-DECAY WITH TACHYONIC 
ANTINEUTRINO 

The tritium /3-decay with tachyonic antineutrino was 
discussed in the papers by Ciborowski and Rembielinski 
pBlp^t . Here for a completeness we present only a brief 
view on results obtained in jl6) . As in the standard case, 
the two graphs shown on the Fig. |l2| contribute to this 
process on the tree level, where dominant is the first one. 
Effectively, for energies much less than mw, it reduces to 



n p + n p + 




FIG. 12. Graphs contributing to the /3-decay 
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the four-fermion interaction with the amplitude square 
given by chiral coupling 



|M ch | 2 = 2G|Tr u^yl-JL^tsy } 
x Tr [upUpj^l - gAl 5 )u n u n j„(l - gAl b )} (54) 

where the axial coupling gA ~ 1.25; in the case of helicity 
coupling, factors (1 — 7 5 )/2 are replaced by 1. The ex- 
plicit calculation of the electron spectrum for chiral and 
helicity couplings as well as comparision with the experi- 
mental data are given in |ll| . Let us stress only that the 
mysterious bump near the end point is reproduced on 
the theoretical plot. More exhaustive discussion of this 
question, as well as of the second anomaly problem in the 
tritium /3-decay can be solved with help of the tachyonic 
hypothesis p6|. 
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IX. CONCLUSIONS 

In this paper we have analyzed the introductory re- 
sults related to tachyonic neutrino hypothesis. We do 
not pretend to explain all the problems with neutrino 
physics by means of this hypothesis. However, from the 
behaviour of the calculated decay rates for three dom- 
inant processes: emission of a vv pair by a neutrino, 
radiative decay v — ► ^7 and /3-decay, it follows that this 
possibility cannot be ruled out of considerations. More- 
over such a hypothesis provides us with an additional 
mechanism simulating flavor oscillations. The calcula- 
tions involving both oscillations and three body decay 
for solar, atmospheric and SN1987A neutrino fluxes are 
now in progress. 
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APPENDIX A: VERTICES IN OUR GAUGE 
MODEL 

The vertices necessary for calculation of amplitudes for 
three processes considered here in our gauge model are 
shown on the Table [II. The rest of the Feynman rules 
in the convention adopted in this paper is given e.g. in 
the Ref. M. 



APPENDIX B: AMPLITUDES FOR NEUTRINO 
DECAY vt -v v t vi,v v FOR ALL COUPLINGS 

The square of amplitude for the chiral coupling 



where 



1 



|M| 2 = 2G%m> iV n lxv 



(uk) 



(Bl) 



1"" = - I 1 , 

2 I \/(uk) 2 + k 2 



(up) 



y/ (up) 2 + K 2 ^ 

(Wp v + k v p^ - (kp)g^ - ie^kapp) 



L (uk) \ 
VM 2 + k2 V V( uk ) 2 + k2 ) 

+ k v u» - (uk)g^ - ie^ afs k a uf3) 

« 2 L , M \ 

2yJ(uk) 2 +K 2 \ x/(u P ) 2 + K 2 J 

+ p v u^ - (up)g^ + ie^ a/3 p a u p ) 



(k»u v 



{p^u v 



and 



2 x /(uk) 2 + ^y/iup) 2 + K' 

v 2 l . funxn ,,2 



(2u^u v - g» v ) 



(ur) 



x/(ul) 2 + fi 2 J \ x/(ur) 2 + y 2 f 
(l^u + - {lr)g^ + is^vaprr 13 ) 

m 2 A , 1 



2 x /(ur) 2 +fi 2 \ VPfVj 
+ IvUn - (uVjg^u + i£,j, va f3l a vP) 

M 2 L + (ur) \ 

2^(ul) 2 + fi 2 \ ^(ur) 2 + t i 2 ) 

+ r u Ufj, - (ur)g^ ~ ie^vapr"^) 



(l»u v 



(r^ 



2y/(ul) 2 + y 2 y/(ur) 2 +n : 



(2u tl u v - g^) 



Squares of amplitude for helicity (s = 0) and 7 s (s = 1) 
coupling read 



|M| 2 = 2G 2 F m^n sfll/ , 



(B2) 



where 



m? = [((-1) s k 2 - kp) 

1 



((-i)V«m 



y (uk) 2 + k 2 y/ (up) 2 + k 2 

-(uk)( U p)) -K A -K 2 ((up) 2 + (Uk) 2 ) 

-(uk)(up)(kp))} 
2k 2 (k 2 - (-l) s (kp)) 

\J (uk) 2 + K 2 y/ (up) 2 + K 2 

(uk)(u P ) - (-l) S K 2 



+ (k» P V + k V P ») 



yj (uk) 2 + K 2 y/ (up) 2 + K 2 



y (uk) 2 + K 2 yj (up) 2 + k 2 
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TABLE III. Vertices in our gauge model used in this work 



Z 



z 



w~ 



H 



v i v 



G° 



v i v 



G~ 



v i t 



G+ ( G+ 



ley 



igj" 1 - 7 5 
2 cos 6w 2 

igY 1 - 7 5 
V2 2 

ign 



2mw 

S- \(m — ft)/ — irn + ftW 5 ! 

2V2m w 1 1 

-ie(p + g) M 



W+ / W+ 




2 cos 9w 



(gv - 3A7 5 ) 





i w+ 


igj" 1 - 7 5 




i , v 


V2 2 




1 if 

1 


igm 




1 I 

— i — . — 


2mw 




, G° 
i 


gmj 5 




! i 
— i — . — 


2mw 


f 


i G+ 
i 

V 

1 V 

— i — . — 


2^2m w 



; [(r - q)\g^ + (q- p)ugx^ + (p- r)^g v \] 



lemwg^v 



+ {u»p u + «y) 



K 2 ((-iy(uk) + {u P )) 



Ct0 



\/(uk) 2 + K 2 ^J(up) 2 + K 2 

K 2 {(~l) S k a + P a ) - {uk)k a pP 
y/(uk) 2 + K 2 



k 2 ((-l) s p a + k a ) ul 3 + {up)k a pP 

\J {up) 2 + K 2 



and 



+ 



= ^(«)[((-l)V + ir) 
1 



((-1) V ((lr) 



+ + Uylp) 



y/(ul) 2 + H 2 y/{ur) 2 + fi 2 

-{ul)(urj) + M 4 + m 2 ({ul) 2 + (m-) 2 ) 
+(«0(«r)(Ir))] 

2M 2 (-/i 2 -(-l) s (/r)) 
^(ul) 2 + fi 2 ^(ur) 2 + fi 2 

(ul)(ur) - (-l)V 



^{ul) 2 + ii 2 ^/(ur) 2 + \x 2 
V 2 ((-iy(ur)-(ul)) 
^(ul) 2 + Li 2 ^(ur) 2 + /j 2 

,i 2 ((-i)'(ui) - H) 



- 1 



v/(m/) 2 + p 2 ^{ur) 2 +/J, 2 



((-l) s r Q - T) + (w)/-"/' 3 
\f{ur) 2 + /x 2 



A1 2 ((-l) s / Q -r Q ) M ' 3 -(u/)r Q Z' 3 



APPENDIX C: AMPLITUDES OF PROCESSES 
CONTRIBUTING TO THE RADIATIVE DECAY 

Amplitudes of processes contributing to the radiative 
decay described by diagrams on Fig. §] (a)-(f) read 



f>f = v(p)M?v(k) 
where i — a,b, c, d, e, / and: 



(CI) 



d A r l + 7 5 „ 
-7 



x 7„- 



2 7 (2tt) 4 2 
[(p — r )l + m]7 a [(/c — r)j + m] 

J r 2 _ to 2 ^] [(j) — r)2 _ m 2j _ r p _ m 2 
1 - 7 5 



(C2) 



ff 2 e /" d 4 r I+7 5 
[r7 + m 



■T 



[r 2 — m 2 ] [{p — r) 2 — m 2 ^] [(k — r) 2 — ™ 2 



wi 



(2p-k- r) x g% + {2r-k- pfgxp 



(C3) 
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g e 



d 4 r 



8m 2 w J (2tt) 



j [(m — k)/ + (m + k)7 5 



Hp ~ r )7 + ™]7 A '[(fc — r)7 + to] 
[r 2 — to 2 ^] [(p — r) 2 — m 2 ] [(fc — r) 2 — to 2 ] 

x [(m - k)I - (m + k)j 5 ] , (C4) 

M '<'£klw? [(m -* ),+(m+ * )li] 



[(p — r)7 + to] (fc + p — 2r)^ [(fc — r)7 + to] 
[r 2 — to 2 ] [(p — r) 2 — to 2 ^] [(/c — r) 2 — to 2 ^] 
x [(to — n)I — (to + «)7 5 ] , 

8 7 (2tt) 4 7 1 7 J 
[r 7 + to] 

[r 2 — TO 2 ][(p — r) 2 — to 2 v ][(/c — r) 2 — to 2 ^] 
x [(to — k)I — (to + k)7 5 ] , 

d 4 



(C5) 



(C6) 



Mi 1 = — / -— t 

' s y (2^)4 



[(to - «)J + (m + k)7 5 ] 
[7*7 + to] 



[r 2 — to 2 ] [(p — r) 2 — to 2 ^] [(fe — r) 2 



x(l+ 7 5 ) 7 r (C7) 

After calculations and renormalization, the integrals con- 
tributing to the form factor Gy in the Eq. (|45|) read 

rl x(x 2 - 3.x + 2) 



/„ = -2 



dx 



Zm 2 w ' 



X 2 K 2 + x{m 2 v — K 2 — TO 2 ) — TO 2 ^ 



x 2 (2x + 1) 



q X 2 K 2 + x(m 2 — K 2 — TO 2 ^) — TO 2 

7 

6771^ ' 



1 X 2 [x(k 2 — TO 2 ) — (TO 2 + K 2 )] 

dx 



L w Jo 



2m 2 



■In 



m 



h = 



rn 



w ->o 



X 2 K 2 + xim 2 ^ — Hi 2 — TO 2 ) — TO 

f in w \ 2 

V TO / 

x(ir, — 1)[x(k 2 — to 2 ) + 2to 2 ] 



2 

ir 



X 2 K 2 + x(m 2 — K 2 — TO 2 ^) — TO 



5to 2 + ft 2 
6to^ 



x 2 k 2 + :e(to 2 — k 2 — mf^) — to 2 



2to 2 k 



APPENDIX D: SUM OVER PHOTON 
POLARIZATION VECTORS 

Sum over photon polarization vectors (on the mass 
shell q 2 = 0) in the CT synchronization reads: 

n""(g)= E eV* 

A=-l 
1 



[(ug)V» 



(wq) 2 

-H(tiY + «V)] 



(Dl) 



Notice, that c^IF" = = 0, so e\, and u^ 1 form 

an pseudo-orthogonal basis. 
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